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Abstract 

Using the mean field theory, a comparative study of the wetting and layering 
transitions of a spin- 1/2 Ising model with perfect and corrugated surfaces, is established. 
The phase diagrams are investigated and compared in the presence of both a longitudinal 
and surface fields. The effect of increasing the temperature and the transverse field on 
the wetting and layering transitions is outlined. 
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1 Introduction 



Recently, the wetting and layering transitions of magnetic Ising systems have been subject 
of great interest. Experimental studies have motivated much theoretical works in order to 
understand and explain the growth of thin layers from only single atoms. A simple lattice-gas 
model with layering transitions has been introduced and studied in the mean field approximation 
by de Oliveira and Griffiths [1] . Multilayer films adsorbed on attractive substrates may exhibit 
a variety of possible phase transitions, as has been reviewed by Pandit et al. [2], Nightingale 
et al. [3], Patrykiejew et al. [4] and Ebner et al. [5-8]. One type of transitions is the layering 
transitions, in which the thickness of a solid film increases discontinuously by one layer as the 
pressure is increased. Such transitions have been observed in a variety of systems including for 
example 4 He [9,10] and ethylene [11,12] adsorbed on graphite. Ebner and Saam [13] carried out 
Monte Carlo simulations of such a lattice gas model. Huse [14] applied renormalization group 
technique to this model. It allowed the study of the effects on an atomic scale in the adsorbed 
layers. The lattice gas models applied to the wetting phenomena was reviewed by Dietrich [15]. 
The effect of finite size on such transitions has been studied, in a thin film confined between 
parallel planes or walls, by Nakanishi and Fisher [16] using mean field theory. 
The model of transverse field was originally introduced by de Gennes [17] for hydrogen-bonded 
ferroelectrics such as KH 2 POa. Since then, this model has been applied to several physical 
systems, like DyV0 2 , and studied by a variety of sophisticated techniques [18-21]. In order to 
study the amorphization of a crystalline Ising ferromagnet with a transverse field, Kaneyoshi [22] 
introduced the differential operator technique as a generalised but approximate Callen relation 
derived by Sa Barreto and Fittipaldi [23]. The system has a finite transition temperature, 
which can be decreased by increasing the transverse field to a critical value fi c . The effect of a 
transverse field on the critical behaviour and the magnetisation curves was studied [18-21] and 
by Kaneyoshi et al. [24,25]. Using the perturbative theory, Harris et al. [26] have studied the 
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layering transitions at T = in the presence of a transverse field. Benyoussef and Ez-Zahraouy 
have studied the layering transitions of Ising model thin films using a real space renormalization 
group [27], and transfer matrix [28] methods. In a previous work [29] we have established the 
effect of temperature and transverse field, for a fixed surface field values, on the wetting and 
layering transitions of an Ising system with perfect surfaces. The purpose of this work is to 
study the effect of a corrugated surface on the wetting and layering transitions, of an Ising 
system spin-1/2, using the mean field theory. The outline of this work is as follows. In Section 
2, the model and method are presented. Section 3 is devoted to study the ground state for 
the perfect and corrugated surface cases. The results and the phase diagrams obtained are 
discussed in Section 4. 

2 Model and method 

We consider a film with finite thickness of a three dimensional spin-1/2 model, with N layers, 
limited by two surfaces. Figure 1 represents the geometry of this system for 

a) two perfect surfaces, 

b) presence of a corrugated surface. 

The Hamiltonian describing this system can be written as 

n = - E Jus's* - + nsf) (i) 

<i,j> i 

where the first summation is carried out over nearest-neighbour pairs of spins. S" = ±1, 
(a = z,x) are the Pauli matrices for a spin i. The exchange interaction = J is assumed to 
be constant, and f2 is the transverse field. The total longitudinal field Hi, applied on a site i, 
is distributed according to the nature of the top surface as it is discussed below. 
For a system with perfect surfaces, Fig. la, Hi is assumed to be uniform in a layer k, and 
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defined by: 



H h = I 



H + H al for k = 1 

H for 1 < k < N 



H + H SN for k = N 



(2) 



where the surface fields H sl and H SN = —H sl are applied on the first layer k — 1 and the last 
layer k = N, respectively. H is the external field applied on the whole system. 
For a corrugated surface film with n steps, each step containing L spins in the x direction, Hi 
is given by: 

H + H s \ for ieSi 
ih = < H for ieS k : 1 < k < N (3) 

H + H SN for ieS N 

where the surface field H sl acts on the corrugated surface Si, and H SN = —H s i is applied on 
the bottom surface S N , Fig. lb. The configurations S k (0 < k < N) are defined in Section 3 for 
both the perfect and corrugated surfaces. H is the external field applied on the whole system. 
The part of the Hamiltonian describing a site i, for the perfect and corrugated surfaces, is given 
by 



n i = (£j ij sf + H i )st + nst 



(4) 



the summation runs over nearest neighbour sites j of %. The diagonalization of the operator Hi 
leads to the eigen values Xf = ±Jxj + Q 2 , where Xi = Y^j^i -hjSf + 



The magnetisation per site can be written as 



rrii = 



TrSfexp(-pHi) Ah + X 



+ \2 



n 2 



Trexp(-pHi) \h + A+) 2 + Q 2 



tanh^xf + n 2 ). 



(5) 



Using the mean field theory, the free energy of the system can be written as follows 

F = ~\ E log(2 cosh(/3A+)) + i E E Jammy 
P i z i j 



(6) 
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For a simple cubic lattice, the magnetisation and the free energy for each plane k are given, 
respectively, by the following expressions: 

m k = ^-[4m k + (m k+1 + m fc _i + H')\ tanh(/3'A fe ) (7) 

and 

F k = ~ log(2 cosh(/3A fc )) + ^m fe (4m fc + (m k +i + m k -i) (8) 

with 

A fe = y/ (4m fc + (m fc+ i + m fe _i) + JT ) 2 + ^' 2 , (9) 
assuming that Jy = J, the reduced parameters are defined by: — f3J, if' — y, Q' — j. 

3 Ground state 

3.1 Perfect surface 

The ground state of this model is established in Fig. 2a. The start point is a situation where 
all spins are down: configuration So. For H si / J < 1, we have only one transition, So <-> Sn- 
By contrast, for H sl / J > 1, we have three transitions namely; the first transition, S <->• Si, 
arises for H/ J = 1 — H sl / J , the second transition, Si <-> Sat-i, occurs for if/ J = 0, while the 
last transition Sjv-i ^ Sn, is seen for if/ J = 1 + i? s i/ ^- The notation Sfc(A; = 0, 1, N) is a 
situation where, starting from the top surface, k layers are spin-up and the remainder N — k 
layers are spin-down. The Fig. 2a shows that the transitions arise for values of H sl j .1 above a 
critical value, H c sl / J = 1. 

3.2 Corrugated surface 

The model is described in section 2, the corresponding ground state is illustrated in Fig. 2b 
for a system formed with iV = 20 layers, n = 2 steps and L — 10 spins, for each step, in the 
x direction. The start point is a situation, called So, where all the spins are down. Under the 
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effect of increasing the surface field H sl /J, and for a bulk field value H/J close to f — H sl /J, 
the spins of the top external surface, called Si, will flip and become up. Increasing the bulk 
field, the spins of the layer 2 will flip, this situation is noted S 2 . Increasing the bulk field more 
and more, the spins of the layer 3 will flip, this configuration is noted S3, and so on. As it is 
shown Fig. 2b, the transition So *-> Si, which appears for H/ J = 1 — H si / J and the transition 
Sjv-i <-> Sat, arising for H/J = — 1 + H s i/J, are fixed only by the surface field H s i/ J value, 
and do not depend on the parameters n and L. By contrast, the number of transitions parallel 
to the line H/ J = is close to the number of corrugated surfaces n, and the corresponding 
value of H/ J is fixed by the number of spins per step L. A generalisation of the ground state 
energy E(Sk, n, L, N), of this system, is given in Appendix A, as a function of arbitrary values 
of the parameters: n, L and N. 

4 Results and discussion 

Numerical results are giving for a system formed with N = 20 ferromagnetic layers of a 
spin— 1/2 Ising model with free bound conditions. In order to examine the temperature effect, in 
absence of transverse magnetic field, fl/ J = 0, on the wetting and layering transitions, we plot in 
Fig. 3 the corresponding phase diagrams. By increasing the temperature we show the existence 
of a critical value, T w / J, above which a sequence of layering transitions occurs this means that 
the temperature effect is sufficient to produce layering transitions. This temperature, T w /J, is 
called the "wetting temperature", which depends on the transverse field and surface magnetic 
field values. The first and last layer transitions occur in the absence of the transverse magnetic 
field, but only under the increasing bulk magnetic field for a temperature T/J> T w /J, the 
other layering transitions arise for increasing temperatures. On the other hand, the corrugation 
effect is to accelerate the layering transitions phenomenon since the transitions S <-> S 2 and 
So <-> S3 are present in the corrugated surface case: Fig. 3b and are absent in the perfect 



6 



surface situation: Fig 3a. The same behaviour is found when inverting the temperature and 
the transverse field roles, as it is summarised in Fig. 4. The phase diagrams show that the 
layer transitions are found for Q/J> Q w /J with increasing transverse magnetic field at fixed 
temperature. fl w /J is called the "wetting" transverse field. The same transitions, found in 
Fig. 3, are still present for the corrugated surface case, Fig. 4b, and disappear for the perfect 
surface situation, Fig 4a. 

Figs. 3 and 4 are plotted for a surface field value H sl /J = 0.9, less than the critical value 1. 
The effect of increasing the surface field above the critical value is summarised in Figs. 5 and 6 
for Hgi/J = 1.2. It is worth to note that the "wetting" parameters, T w j J and Q w /J, disappear 
since the layering transitions are present even for very low temperature: Figs 5a and 5b, and 
very low transverse field: Figs 6a and 6b. Furthermore, the transitions S 1 <-> S 2 and S 2 <-> S 3 
are present in Figs. 5b and 6b for T/J — > and Vt/J — >• 0, respectively, in good agreement 
with the ground state illustrated in Fig. 2b, for the corrugated surface case. On the other 
hand, the dependency of the wetting temperature T m j J on the wetting transverse field VL W / J, 
is illustrated in Fig. 7, corresponding to the surface field value H s i/J = 0.9, for the both cases: 
perfect and corrugated surfaces. As it has been outlined above, the corrugation effect is to 
accelerate the wetting phenomenon since the wetting region, of the phase diagram in Fig. 7, 
occurs first for the corrugated surface before that one corresponding to the perfect surface. 

5 Conclusion 

Within the mean field theory, we have studied the wetting and layer transitions of a 3-D 
spin— 1/2 Ising model. For a surface field less than the critical value 1, the system exhibits 
a sequence of layering transitions only above a wetting temperature (wetting transverse field) 
for fixed transverse field (temperature). For surface fields higher than the critical value, the 
layering transitions under only under the effect of increasing the bulk field even for Q/ J = or 
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T/J^O. Furthermore the corrugation effect is accelerate the layering transitions and wetting 
phenomena compared to perfect surfaces. 
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Figure Captions 

Figure 1. Geometry of the model with N layers in the z direction. The system is infi- 
nite in the direction y with: 

a) perfect surfaces case, 

b) presence of a corrugated surface with n steps, each step contains L spins in the x direction. 
Figure 2. Ground state phase diagrams in the (H/ J, H sl / J) plane for: 

a) perfect surfaces case with N = 20 layers, 

b) corrugated surface case with N = 20 layers, n = 2 steps and L — 10 spins for each step. 
Figure 3 Phase diagrams in (H/J,T/J) plane for H s i/J = 0.9 and Q/J = 

a) perfect surfaces case, 

b) presence of a corrugated surface. 

Figure 4. Phase diagrams in (H/J,Q/J) plane for H s i/J = 0.9 and T/J — 0.5 

a) perfect surfaces case, 

b) presence of a corrugated surface. 

Figure 5. Phase diagrams in (H/J,T/J) plane for H sl / J = 1.2 and Q/J = 

a) perfect surfaces case, 

b) presence of a corrugated surface. 

Figure 6. Phase diagrams in (H/ J,Vt/ J) plane for H s i/J = 1.2 and T/J = 0.5 

a) perfect surfaces case, 

b) presence of a corrugated surface. 

Figure 7. The wetting temperature T w / J dependency on the wetting transverse field Q w / J, 
with fixed surface field H s \/ J = 0.9, for the both cases: corrugated and perfect surfaces. 
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Appendix 

Expression of the energy E(Sk, n, L, N): 

In the following we will investigate, the ground state energy E(Sk, n, L, N) as a function of the 
parameters S k , n, L and N (see text). For T = the Hamiltonian of Eq. (1) allows to write : 

E(S k , n, L, N)/J = A x H sl /J + Bx H/J + C. (10) 

Where A = N Hsl is the coefficient of H sl , B = N H is the coefficient of H and C = Nj is the 
coefficient of J, which depend on the parameters Sk,n,L and N 

For a given configuration Sk,(x = 0,1, N), of a corrugated surface system formed with 
N layers, n steps with L spins for each step. These expressions will generate a sequence of 
transition lines from a configuration Sk to a configuration Si when < k, I < N. 

.1 Coefficients Nn sl (Sk, n, L, N) 

1- Configuration So (all spins down): NH 3l (So,n, L, N) = 2n. 

2- Configuration S k , (k = 1, N - 1)): N Hal (S k ,n, L, N) = -2(n + 1)(2L - 1). 

3- Configuration Sn ( all spins up): Nn sl (SN,n, L, N) = —In. 

.2 Coefficients N H (S k , n, L, N) 

1- Configuration So: 

N H (S ,n, L, N) = YT z i{{2z{L - 1) + 1) + (N - (n + l))(2(n + 1)(L - 1) + 1). 

2- Configuration Sj : 1 < j < n + 1: 

7Vtf(^, n, L,N) — — El=i(2z(L - 1) + 1) + £S+i(2*(L - 1) + 1 - 2L) + (iV - (n + l))(2(n + 
1)(L-1) + 1). 

3- Configuration Sj : n + 1 < j < N: 

N H (Sj, n, L,N) = - J2 n z tl(2z(L - 1) + 1) - (j - (n + l))(2(n + 1)(L - 1) + 1) + (N - j)(2(n + 
1)(L-1) + 1). 
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.3 Coefficients Nj(Sk, n, L, N) 

1- Configurations S : 

Nj(S ,n,L,N) = -(E^(4(2z(L - 1)) + 3) + (N - 1 - (n + l))(4(2(n + 1)(L - 1)) + 3) + 
(3(2(n+l)(L-l)) + 2)). 

2- Configuration S^: 

Wj^i, ra, L, N) = -((4L + 1) + £? =2 (2(4 + 2(L - 2) + 4(z(L - 1) - L)) + 3) + (2(2(L - 1) + 
4((n + 1)(L - 1) - L)) + 3) + Ef=„ 1 +2 (4(2(n + 1)(L - 1)) + 3) + (3(2(n + 1)(L - 1)) + 2)). 

3- Configuration Sj : 1 < j < n + 1: 

JVj^-, n, L, N) = -(Et 1 1 (5(2z(L - 1) + 1) - 1) + E^(2(4 + 2(L - 2) + 4(z(L - 1) - L)) + 
3) + (2(4 + 2(L - 2) + 4((n + 1)(L - 1) - L)) + 3) + (2(2(L - 1) + 4((n + 1)(L - 1) - L)) + 3) + 
Ef=- 1 +2 (4(2(n + 1)(L - 1)) + 3) + (3(2(n + 1)(L - 1)) + 2)). 

4- Configuration S' n+ i: 

Nj(S n+1 ,n, L, N) = -(ELi(5(2^(L - 1) + 1) - 1) + (2(2(n+ 1)(L - 1)) + 1) + Ef =n 1 +2 (4(2(7i + 
l)(L-l)) + 3) + (3(2(n + l)(L-l)) + 2)). 

5- Configuration £j : n + 1 < j < N — 1: 

Nj(S^n,L,N) = -(E: =1 (5(2z(L - 1) + 1) - 1) + E^i(4(2(n + 1)(L - 1)) + 3) + (2(2(n + 
1)(L - 1)) + 1) + Ef = ~ 1 + i(4(2(n + 1)(L - 1)) + 3) + (3(2(n + 1)(L - 1)) + 2)). 

6- Configuration 5^: 

Nj(S N ,n,L,N) = -(£? =1 (5(2z(L - 1) + 1) - 1) + E^~i 1 (4(2(n + 1)(L - 1)) + 3) + (3(2(n + 
l)(L-l)) + 2)) 
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